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Abstract 

Wc report preliminary results of AKARI observations of two globular clusters, NGC104 and NGC362. 
Imaging data covering areas of about 10 x 10 arcmin 2 centered on the two clusters have been obtained 
with InfraRed Camera (IRC) at 2.4, 3.2, 4.1, 7.0, 9.0, 11.0, 15.0, 18.0 and 24.0 /xm. We used F u /F 2 and 
F24/F7 flux ratios as diagnostics of circumstellar dust emission. Dust emissions are mainly detected from 
variable stars obviously on the asymptotic giant branch, but some variable stars that reside below the tip 
of the first-ascending giant branch also show dust emissions. We found eight red sources with F24/F7 ratio 
greater than unity in NGC362. Six out of the eight have no 2MASS counterparts. However, we found no 
such source in NGC104. 

Key words: infrared:stars stars:AGB and post-AGB Galaxy: globular clusters individual (NGC104, 
NGC362) 



1. Introduction 

In order to calibrate the geometric distortion of the 
InfraRed Camera (IRC; Onaka et al. 2007) onboard 
AKARI (Murakami ct al. 2007), we observed two galac- 
tic globular clusters, namely NGC104 (= 47 Tuc) and 
NGC362 during the AKARI in-orbit performance veri- 
fication phase. Four pointings were dedicated for these 
observations. Here, we will use these data to study the 
evolution of low-mass stars with an emphasis on mass loss 
from them. 

Galactic globular clusters are the best test grounds for 
stellar evolution theories of low-mass stars, since we can 
reasonably assume that they are composed of a single 
stellar population - that is, their constituent stars were 
formed at the same time, in the same volume of space, 
and from the same cloud of gas. Understanding of the 



evolution of low-mass stars is important because they may 
lose as much as ~ 40% of their initial mass during their 
life time (Wachter et al. 2002), and also, as they rep- 
resent the majority of the stars in the Galaxy, their role 
in the galactic formation and evolution is not negligible 
(e.g., Schroder and Scdlmayr 2001). 

The mass loss process, which dominates the evolution 
of the star itself, is still a poorly understood phenomenon. 
Many observations (e.g., Tanabe et al. 1997) showed that 
the phase of the highest mass loss rate is achieved dur- 
ing the last stage of the asymptotic giant branch (AGB) 
evolutionary phase. Although it is not a large scale, how- 
ever, mass loss also occurs along the first-ascending giant 
branch (RGB). It should exert a greater effect on stellar 
evolution, especially for low mass stars, given that their 
residence time in the RGB is much longer than that in the 
AGB (Schroder and Cuntz 2005). 
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NGC104 and NGC362 are well suited for studying a 
mass loss history along the RGB and AGB, because of 
well populated RGB and AGB stars. AKARI/IRC can 
detect red giants well below the tip of the RGB in the two 
clusters. The basic parameters of NGC104 and NGC362 
are listed in Table 1, and a short description of each cluster 
including introductions of previous work follows. 

Table 1. The basic parameters of NGC104 and NGC362. 





NGC104 


NGC362 


DM [mag]* 


13.50 ± 0.08 1 


15.06 3 


[Fe/H] 


-0.66 2 


-1.12 3 


Age [Gyr] 


11.2 ± l.l 1 


8 ~ 9 3 



* Distance modulus. 

References: 1 Gratton et al. (2003), 2 Carretta and 
Gratton (1997), 3 Gratton et al. (1997) 

1.1. NGC104 

NGC104 is the second brightest (after u> Cen) galac- 
tic globular cluster. To date, 42 long period variables 
have been found (Lebzelter & Wood 2005). Mid-IR sur- 
veys of this cluster have been done by using the ISOCAM 
(Cesarsky et al. 1996) on board the Infrared Space 
Observatory (ISO; Kessler et al. 1996) to study the 
evolution of dust mass loss along the RGB and AGB 
(e.g., Ramdani and Jorissen 2001, Origlia et al. 2002). 
Ramdani and Jorissen (2001) observed outer regions of 
NGC104, and Origlia et al. (2002) observed an area of 
about 5' x 3' centered on the cluster core. The IRC's 
wide field of view (10' x 10') enables us to cover both the 
central and outer regions within a single observation. 

1.2. NGC362 

This cluster is also well studied. Along with NGC288 it 
forms one of the most famous "second parameter" couples. 
NGC362 and NGC288 have about the same metallicities 
(= "the first or main parameter"), but their horizontal- 
branch (HB) morphologies are different. NGC362 has a 
red HB morphology, while NGC288 has a blue one (e.g., 
Sandage and Wallerstein 1960, Bolte 1989, Green and 
Norris 1990). To explain the difference, there must be 
at least one " second parameter" . There are many sec- 
ond parameter candidates, e.g., cluster age, mass loss 
along the RGB, helium abundance, rotation and deep 
helium mixing, dynamical interactions involving binaries 
and even planets, environmental effects in high-density 
environments, and so on (e.g., Vink and Cassisi 2002 
and references therein), and its origin is still a controversy 
(e.g., Stetson et al. 1996). Mass loss episodes along the 
RGB may be an "inclusive" second parameter (Catelan 
et al. 2001), as they are derivative (i.e., determined by 
other parameters of the star). Origlia et al. (2002) also 
observed this cluster, finding three stars with mid-infrared 
excess. 



2. Observations & Reductions 

Imaging observations of NGC104 and NGC362 were 
obtained on 2006 May 1st (NGC104), and 6th and 7th 
(NGC362) UT with the IRC. The IRC03 AKARI IRC 
observing template (AOT03) was used, yielding imaging 
data at 2.4, 3.2, 4.1, 7.0, 9.0, 11.0, 15.0, 18.0 and 24.0 /xm 
taken in at least 2 dithered positions. 

Raw data were processed with the IRC imaging data 
pipeline, version 070104 (see IRC Data User's Manual 
Lorentc et al. (2007) for details). The resultant IRC mo- 
saic images have pixel sizes of l."46 pixel -1 , 2."40 pixel -1 
and 2"38 pixel" 1 for NIR (2.4, 3.2, and 4.1 /mi), MIR-S 
(7.0, 9.0, and 11.0 /mi) and MIR-L (15.0, 18.0 and 24.0 
/im) channels of IRC, respectively, covering an area ~ 100 
arcmin 2 around the cores of NGC104 and NGC362 in each 
wavelength. 

2.1. Photometry 

To derive calibrated fluxes for each star, point spread 
function (PSF) fitting photometry was performed on the 
mosaiced images with the IRAF 1 package DAOPHOT. 
Photometry was done for each mosaiced image indepen- 
dently. This involved the following steps: 

1. DAOFIND was used to find stars whose fluxes are 
at least 5 a above the background, where a is the 
background noise estimated locally around stars. 

2. Aperture photometry was performed on all of the 
stars found in step 1, using the task PHOT with 
aperture radii of 10.0 and 7.5 pixels for NIR and 
MIR— S/MIR— L images, respectively. We used the 
same aperture radii as had been used in the standard 
star flux calibration (Tanabe et al. in preparation), 
so the aperture corrections were not applied. The 
resultant astronomical data units were converted to 
the calibrated fluxes by using the IRC flux calibra- 
tion constants version 070119. 

3. Several stars with moderate flux (i.e., with a good 
signal-to-noise ratio and unsaturated) and without 
neighbors within 7 pixels were selected from the re- 
sults of step 2. We found more than 5 such stars in 
each mosaiced image. The selected stars were used 
to construct a model PSF. 

4. The PSF fit was adopted to all of the stars found 
in the mosaiced images using ALLSTAR to get 
their instrumental fluxes and their corresponding er- 
rors. To check the array-location-dependence, we 
ran ALLSTAR with an option that the PSF can be 
linearly variable over the images. With this test, we 
found that the PSF does not vary significantly over 
the array. Hence constant PSF is assumed over an 
image. 

5. The resultant instrumental fluxes were shifted so 
that the instrumental fluxes of the stars selected in 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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step 3 match the calibrated fluxes calculated in step 
2. 

We do not deredden the measured fluxes, but we applied 
color corrections on the calibrated fluxes by assuming a 
black body with the effective temperature of 3500 K. The 
temperature change of ± 500 K yields 2.4% change in the 
correction factor in the 2.4 /xm case, but less than 1.0% 
for the others. Therefore, any discussions followed are 
almost insensitive to the assumed reference black body 
temperature. 

2.2. Cross-identification with 2MASS sources 




* * 



ric coordinates of detected sources in the mosaiced images 
of each wavelength. As the result, their coordinates were 
determined to the accuracy of about 2"0 relative to the 
corresponding 2MASS sources. In Figure 1, we show the 
results of the cross-correlation between AKARI IRC 11.0 
jum and 2MASS sources. The background images are IRC 
11.0 fim images of NGC104 (top) and NGC362 (bottom). 
Crosses represent the IRC 11.0 /jm sources with 2MASS 
counterparts within a radius of l."0. 

After this cross-correlation, we inspected each source 
with -fr'2MASS < 13.0 mag by the eye to remove miss- 
identifications. Since we worked with the photometric 
data in high stellar density regions, we eliminated any 
IRC sources that were merged (we assumed that any IRC 
sources that have neighbor(s) within a half width at half 
maximum of the PSF of each wavelength as merged) with 
a brighter star, because they would introduce elements of 
confusion into the following discussions. 
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Fig. 1. AKARI IRC 11.0 /im image of NGC104 (top) and 
NGC362 (bottom). Crosses show IRC 11.0 (im sources with 
2MASS counterpart. North is up, and east is to the left. 

Cross-correlation of the IRC sources with 2MASS 
sources has been made in order to determine the astromet- 



Fig. 2. The F11/F2 vs K B diagram for sample stars with 
the variable stars identified. We use the nomenclature of 
Lebzelter & Wood (2005) for NGC104, and Szekely et al. 
(2007) for NGC362, except for two stars C2 and C16 (see 
text). The F2 flux densities were calculated from K B magni- 
tude by adopting the zero-magnitude flux of 666.7 Jy (Cohen 
et al. 2003). See section 3.3 for the two triangles in NGC362. 
Closcup around the crowded part is shown in the inset. 
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3. Results & Discussion 

3.1. Infrared excess of cluster stars 

Ramdani and Jorissen (2001) showed that the ratio of 
the ISO 11.5 /mi to DENIS K s 2 /Am flux density is a 
good indicator of dust mass loss. Therefore, we made a 
IRC 11.0 /un to 2MASS K s flux density ratio (Fu/F 2 ) 
vs ^2MASS diagram of our sample to see which stars 
show circumstellar dust emission. Hereafter, we denote 
2MASS K s as K s unless otherwise described. The top 
panel of Figure 2 is plotted for NGC104, and the bottom 
for NGC362. The K s fluxes (F 2 ) were calculated from K s 
magnitudes by adopting the zero-magnitude flux of 666.7 
Jy (Cohen et al. 2003). We did not apply reddening cor- 
rection on 2MASS nor IRC fluxes. To get a basic idea, we 
calculated the Fn/F 2 ratio expected for a dust-free stellar 
atmosphere by using ATLAS9 (Kurucz 1993; Sbordone 
et al. 2004; Sbordone 2005). It came out that a star 
of T e fj = 3500 K, log g = 1.5, = 2.0 km/s, and 

[Fe/H] = -1.0 would have F u /F 2 ratio of about 0.065, 
showing that the Fn/F 2 ratios for most of the stars in 
NGC104 and NGC362 are consistent with the value pre- 
dicted for dust-free photospheres. Further experiments 
with ATLAS9 showed that metallicity has a negligible im- 
pact on the Fn/F 2 ratio. 

As it has been suggested by several authors (e.g., 
Vassiliadis and Wood 1993), stellar pulsations play a 
key role to trigger dust mass loss. The present results 
also show that all of the bright red giants with Fn excess 
(F11/F2 > 0.1) are indeed variables. Throughout this pa- 
per, we use the nomenclature of Lebzelter & Wood (2005) 
for variable stars in NGC104, and of Szekcly et al. (2007) 
for ones in NGC362, except for C2 and C16. We have to 
note that the variable stars C2 and C16 in NGC362 are 
not listed in Szekely et al. (2007), but they are identical 
to V2 and V16 listed in Clement et al. (2001). It is con- 
firmed by a near-infrared monitoring survey conducted by 
one of us (N. Matsunaga) that C2 is a semi-regular vari- 
able with a pulsation period of about 90 days, and also 
that C16 is a Mira-like variable with a pulsation period of 
about 135 days (Sawyer 1931; Matsunaga 2007). 

Two sources in NGC104 with K s ~ 12.0 and one source 
in NGC362 with K s ~ 13.4 may have Fn/F 2 ratios greater 
than 0.1, but their photometric errors prevent us from re- 
garding them as sources with infrared excess. Flux deter- 
minations of these three stars need further examination 
because they have much larger error bars than other stars 
at similar magnitudes. Put it all together, we found eleven 
stars (VI, V2, V3, V4, V8, V13, V21, LW7, LW10, LW11, 
and LW13) and two stars (C2 and C16) with infrared ex- 
cess (F11/F2 > 0.1) in our observed fields of NGC104 and 
NGC362, respectively. 

3.2. Infrared excess from RGB stars in NGC104? 

Figure 2 also illustrates that a few variable stars below 
the tip of the first red giant branch (TRGB) in NGC104 
do exhibit infrared excess. The TRGB occurs at K = 6.75 
- 7.1 mag in NGC104 (Fcrraro et al. 2000; Lebzelter & 
Wood 2005), and specifically, V13, LW7, and LW11 show 



large (Fn/F 2 > 0.1) infrared excess although they reside 
below the TRGB. These three stars have similar pulsation 
periods of about 40 days and also similar (J — K) colors 
of about 1.0 (Lebzelter & Wood 2005). The period and 
(J — K) color are both being typical of variable stars be- 
low the TRGB found in the Large and Small Magellanic 
Clouds (Kiss and Bedding 2003; Ita et al. 2004a; Ita et 
al. 2004b). As in Ita et al. (2002), stars below the TRGB 
could be either on the AGB or the RGB, but a substantial 
fraction could be RGB stars. The definite identification 
of the evolutionary stage of these stars is difficult based 
solely on the present data. We cannot rule out the pos- 
sibility that they are thermally-pulsating AGB stars dur- 
ing the phase of quiescent helium burning, when stars are 
fainter by just over a magnitude (Marigo et al. 2003). 
Further observations of each star is definitely needed for 
detailed study. 
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Fig. 3. The F24/F7 vs F7 diagram for sample stars with the 
variable stars identified as in figure 2. The dashed lines shows 
5 a sensitivity limit for AOT03 in one pointed observation, 
as given in Onaka et al. 2007. Note that fluxes are color- 
corrected, but not dereddened. Closeup around the crowded 
part is shown in the inset. 

Interestingly, V13 in NGC104 has a F24/F7 ratio com- 
parable to that of VI (see the top panel of Figure 3), which 
has the longest pulsation period, largest pulsation ampli- 
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Fig. 4. The spectral energy distribution of V13 (filled cir- 
cles), 2MASS 00234761-7202498 (asterisks), and VI (tri- 
angles) in NGC104. Note that the fluxes of 2MASS 
00234761-7202498 and VI are scaled so that the H-band 
flux densities of them are matched to that of V13. 
Photometric errors are smaller than the size of the marks. 

tudc, and highest luminosity among all the members of 
NGC104, and shows the typical silicate dust feature in its 
mid-IR spectrum with a mass-loss rate of ~ 10~ 6 M Q yr _1 
(van Loon ct al. 2006). We show the spectral energy dis- 
tribution (SED) of V13 in Figure 4 using 2MASS JHK S 
and all of the IRC fluxes. To emphasize the uniqueness 
of V13 compared to other normal red giants and also to 
genuine mass-losing AGB stars, the SED of a red giant, 
namely 2MASS 00234761-7202498 in NGC104 that has 
similar K s magnitude and (J — K s ) color (K s = 7.606 
mag, and J — K s = 1.081 mag) as those of V13 (K s = 
7.755 mag, and J — K s = 1.089 mag), and also VI are in- 
cluded in the figure. For comparison, the fluxes of 2MASS 
00234761-7202498 and VI arc multiplied by a factor of 
1.022 and 0.410, respectively, that were calculated as the 
ff-band flux densities of V13, 2MASS 00234761-7202498, 
and VI are matched. We also calculated the absolute 
bolomctric magnitude of V13 by fitting two black body 
curves with T eff = 3261 K and T eff = 263 K to the ob- 
served data. We obtained Mu i ~ —3.11 mag as a result, 
using the distance modulus of 13.5 mag for NGC104 (see 
Table 1). Lcbzcltcr ct al. (2006) took low-resolution mid- 
infrared (7.6—21.7 /im) spectra of V13 with the Spitzer 
telescope. They showed that V13 is devoid of a 9.7 /tin 
emission band feature of amorphous silicate, but it has 
broad emission features at 11.5 /im (likely to be AI2O3), 
13 /xm (likely to be an Al— O stretching vibration), and 
20 /im (no firm identification). Our results are consistent 
with theirs, showing that the infrared excess is detected 
at longward of 11.0 /im (Figure 4). Aluminium oxide fea- 
tures have been detected from low mass loss rate oxygen- 
rich AGB stars (Onaka et al. 1989; Kozasa and Sogawa 
1997). Combined with the fact that the Fn/F 2 ratio of 
V13 is not so large, these results show that dust compo- 
sition of V13 is different from those of usual mass losing 
AGB stars. 



3.3. Very red sources in NGC362 

Boyer et al. (2006) observed M15 with IRAC (Fazio et 
al. 2004) and MIPS (Rieke et al. 2004) onboard Spitzer 
Space Telescope (Werner et al. 2004). Their F 2 4/F 8 vs 
Fg diagram revealed that there are at least 23 red sources 
in M15. They suggested that the red sources are mass- 
losing AGB or post- AGB candidates after consideration 
of their loose spatial distribution, 

We use our 7.0 and 24.0 /im data to make the F24/F7 
vs F7 diagram, as shown in Figure 3. The top panel is 
for NGC104 and the bottom is for NGC362. It is seen 
that there are no sources with F24/F7 > 1.0 in NGC104, 
but there are eight red sources in NGC362. Six out of 
the eight have no 2MASS counterparts within a radius 
of 6". The other two sources with 2MASS counterparts 
have large F11/F2 ratios, as seen in Figure 2 (triangles). 
The spatial distribution of the eight sources are shown in 
Figure 5, indicating that they are distributed around the 
cluster, and are not biased to the cluster center. 
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Fig. 5. AKARI IRC 24.0 fim image of NGC362 with eight 
red sources (F24/F7 > 1.0, sec text) identified, their spectral 
energy distributions are shown in Figure 6. North is up, and 
east is to the left. 

According to a model calculation (Pearson et al. in 
preparation), we can expect one or two galaxies in the 
IRC 100 arcmin field of view down to the 5 a sensitiv- 
ity limit of 24 /im. Then, there is little possibility that 
these sources are all background galaxies. Absence of red 
sources in NGC104 also suggests it unlikely that they are 
all galaxies. It is also unlikely that they are mass-losing 
AGB stars in NGC362, because their F7 flux densities 
are too faint (Groenewegen 2006). Boyer et al. (2006) 
found similar objects in M15 by Spitzer observations and 
pointed out the possibility that they could be post- AGB 
stars. However, the fact that no clear counterparts have 
been seen at NIR wavclcgnths may be incompatible with 
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Fig. 6. Spectral energy distributions of the eight red sources 
in NGC362 showing K s (if available) and all of the 
AKARI/IRC data. Note that some stars were not detected 
even in the NIR channel of IRC. The scales on the x- and 
y-axis are the same for each panel. The labels are as in Figure 
5. 

the post-AGB identification. NGC362 lies near the Small 
Magellanic Cloud in projection, therefore they might be 
bright high-mass-losing AGB stars in the SMC. However, 
such stars should be very rare, and we would not expect 
to detect eight of such stars in a 100 arcmin 2 field. We 
show the spectral energy distributions of the eight sources 
in Figure 6. It can be seen that infrared excess is de- 
tected not only at 24.0 /xm but also at 15.0 and 18.0 /im. 
Therefore it is likely that the excess is continuum emis- 
sion and that the contribution from line emission such as 
[O IV] 26 /im is insignificant. Some spectra show peculiar 
features, but the presence of the excess seems to be se- 
cure (based on the multi-band photometry). To identify 
the eight red sources and also to confirm the photometric 
results, infrared spectroscopic observations with AKARI 
and/or Spitzer would be needed. 

4. Summary 

We presented the first AKARI/IRC imaging data of 
NGC104 and NGC362. We found that stars with large 
F11/F2 ratio are all variables, confirming the close link be- 
tween mass loss and stellar pulsation. We detected eight 
sources with large F24/F7 ratio in NGC362, but no such 
sources were present in NGC104. We also showed that 
V13 in NGC104 has interesting features, as being a pos- 



sible RGB candidate with infrared excess. 
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